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PREFACE 


This  report  describes  research  performed  by  S-CUBED  as  part  of  a 
continuing  effort  to  develop  improved  analytical  methods  for  modeling 
the  constitutive  behavior  of  concrete.  The  work  was  jointly  supported 
by  the  Defense  Nuclear  Agency  (Contract  DNA001-84-C-0127)  and  the  Air 
Force  Office  of  Scientific  Research  (Contract  F49620-84-C-0029) .  The 
Principal  Investigator  for  both  contracts  was  Dr.  H.  E.  Read,  and  the 
corresponding  Contract  Technical  Monitors  were  Mr.  C.  B.  McFarland 
(DNA)  and  Lt.  Col.  L.  D.  Hokanson  (AFOSR) . 

Dr.  K.  C.  Valanis,  Professor  of  Theoretical  Mechanics  at  the 
University  of  Cincinnati,  participated  in  the  investigation  through 
his  continuing  role  as  a  consultant  to  S-CUBED. 

The  authors  express  their  appreciation  to  Dr.  D.  H.  Brownell, 
Jr.,  and  Mr.  R.  G.  Herrmann  who  provided  excellent  computational 
support  throughout  the  course  of  the  investigation. 
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SECTION  1 
INTRODUCTION 

The  ability  to  predict  the  behavior  and  strength  of  reinforced 
concrete  structures  is  often  crucial  for  design  purposes  and  depends 
very  strongly  on  our  understanding  of,  and  capability  to  analytically 
model,  the  nonlinear  response  of  plain  concrete  to  applied  stresses. 
From  a  constitutive  modeling  viewpoint,  plain  concrete  is  an 
exceedingly  complex  material  whose  constitutive  properties  under 
general  load  paths  have,  to  date,  largely  eluded  researchers.  As 
noted  by  Chen  and  Suzuki, (*)  a  general  constitutive  relation  capable 
of  describing  both  pre-  and  post-failure  behavior  of  plain  concrete 
simply  does  not  exist  today.  In  fact,  a  general  model  capable  of 
describing  the  behavior  of  plain  concrete  under  all  possible 
conditions,  and  yet  is  sufficiently  simple  for  use  in  design 
applications,  may  be  out  of  reach  at  the  present  time. 

In  recent  years,  a  variety  of  nonlinear  constitutive  models  has 
been  proposed  in  the  literature  for  describing  the  response  of  plain 
concrete  to  short-term  loads.  These  models  have  been  based  upon 
nonlinear  elasticity,  ^)  classical  plasticity,^)  hypo-elasticity,^) 
endochronic  plasticity,^'®)  bounding  surface  plasticity,^'®) 
continuous  damage  theory,  plastic-fracture  theory^^)  and 
microstructura I  considerations.^^)*  In  almost  all  cases,  the  r„odels 
were  developed  from,  and  validated  against,  data  from  standard 
laboratory  tests  in  which  at  least  two  of  the  principal  stresses  were 
equal.  Some  true  triaxial  experiments  have  also  been  performed  on 
concrete,  but  these  have  been  largely  limited  to  proportional, 

*.  For  a  rev i ew  and  critique  of  various  modeling  approaches,  see  Chen 
and  Ting. t13) 


monoton  i  ca  I  I  y- i  ncreas  i  ng  loading  to  f  a  i  I  ure .  the  results  from 

these  tests  are  inadequate  for  identifying  and  characterizing  the 
salient  constitutive  features  of  concrete  under  arbitrary  load  paths, 
unloading,  stress  reversals  and  reloading,  which  are  generally 
encountered  in  practice. 

To  fill  the  need  for  data  on  the  response  of  plain  concrete  to 
complex,  multi-axial  loading  paths,  Gerstle,  Ko  and  coworkers  recently 
conducted  an  extensive  laboratory  testing  program, using  a  true 
triaxial  test  device  designed  specifically  for  use  in  studying  the 
multi-axial  behavior  of  geomaterials.  This  device  independently 
applies  three  normal  stresses  of  up  to  15,000  psi  to  the  sides  of  4- 
inch  cubical  specimens.  The  focus  of  this  experimental  study  was 
confined  to  the  range  of  load  levels  below  that  at  which  significant 
cracking  of  the  plain  concrete  would  occur.  A  total  of  47  different 
non-standard  stress  paths  were  investigated.  Of  these,  a  number  of 
paths  were  chosen  to  provide  insight  into  possible  stress-path 
dependence . 

From  the  results  of  this  study,  it  was  found  that  plain  concrete 
exhibits  behavior  under  complex  loading  which  appears  to  be  not  well 
represented  by  any  of  the  current  models.  In  particular,  it  was  found 
that  the  loading  surface  has  the  form  of  a  box  in  principal  stress 
space,  with  sides  parallel  to  the  principal  stress  axes;  hence,  the 
projection  of  the  loading  surface  on  the  deviatoric  plane  is  a 
triangle,  which  is  aligned  with  the  principal  deviatoric  stress  axes 
in  the  deviatoric  plane  so  that  each  of  its  sides  is  perpend i cu I ar  to 
one  of  the  principal  axes.  Thus,  the  loading  surface  is  not  affine 
with  the  failure  surface.  Except  near  the  failure  surface,  the 
plastic  flow  law  appears  to  be  associated,  with  the  plastic  strain 
rate  vector  normal  to  the  loading  surface  at  the  current  stress  point 
As  the  failure  surface  is  approached,  the  plastic  strain  rate  vector 


begins  to  deviate  from  normality  and  ultimately  becomes  normal  to  the 
fa i I ure  surface . 

In  an  effort  to  analytically  represent  the  results  from  their 
experimental  study,  Stankowski  and  Gerstle  developed  a  simple,  hypo- 
elastic-type  constitutive  model that  is  an  extension  of  a  previous 
model  by  Gerstle. The  model  was  applied  to  several  of  the  complex 
axisymmetric  stress  paths  studied  experimental ly  and,  considering  the 
simplicity  of  the  model,  it  showed  remarkable  predictive  capability. 
Several  other  more  complex  concrete  models^'^®)  were  also  exercised 
around  the  same  paths  by  Gerstle  and  Willam,^^  but  they  failed  to 
show  comparable  predictive  capabilities.  While  the  simple  model  has 
shown  its  ability  to  capture  many  of  the  response  features  of  plain 
concrete  under  several  axisymmetric  stress  histories,  it  has  not,  to 
date,  been  extended  to,  nor  validated  against,  more  general  loading 
conditions.  In  order  to  capture  certain  features  of  behavior  that  are 
beyond  the  scope  of  the  simple  theory,  and  to  have  a  theory  that  is 
applicable  to  general  multi -axial  loading  states,  we  explored  the 
potential  offered  by  the  new  endochronic  theory  for  describing  the 
extensive  experimental  results  given  in  Reference  17. 

The  new  endochronic  plasticity  theory  was  originally  developed 
by  Valanis,^*)  and  has  since  been  applied  with  remarkable  success  to 
various  problems  of  metal  pi ast i c i ty and  geomater i a  Is. The 
theory  is  based  upon  the  hypothesis  that  the  current  state  of  stress 
in  a  material  is  a  linear  functional  of  the  entire  history  of 
deformation,  with  the  history  defined  with  respect  to  a  time  scale, 
called  intrinsic  time,  which  is  itself  a  property  of  the  material  at 
hand.  Such  an  approach  does  not  require  the  notion  of  yield  surface 
nor  the  specification  of  unloading-reloading  criteria,  and  these 
unique  features  make  the  theory  particularly  attractive  for  describing 


the  behavior  of  concrete,  which  does  not  exhibit  a  well-defined  yield 
point.  It  should  be  noted  that  the  new  endochronic  plasticity  theory 
considered  herein  is  substantially  different  from  those  versions  of 
the  older  theory  used  to  develop  the  concrete  models  described  in 
References  5,  6,  and  27.  Proper  closure  of  hysteresis  loops  is 
guaranteed  in  the  present  theory,  so  that  artifices,  such  as  the 
jump-kinematic  hardening  introduced  into  the  concrete  model  described 
in  Reference  27,  are  not  required.  Some  of  the  basic  inelastic 
response  features  of  this  new  theory  have  recently  been  discussed  by 
Trangenstein  and  Read^^)  and  Murakami  and  Read.^®) 


The  purpose  of  the  present  report  is  to  formulate  an  endochronic 
plasticity  theory  for  plain  concrete,  based  upon  the  general  framework 
described  in  Reference  21,  and  to  explore  the  capabilities  of  this  new 
theory  to  predict  the  wide  range  of  experimental  results  for  complex 
multi-axial  load  paths  reported  in  Reference  17.  For  this  purpose, 
the  model  is  first  fit  to  a  small  subset  of  the  data,  and  then 
exercised  around  a  variety  of  different  non-standard  stress  paths  not 
used  in  fitting  the  model.  It  is  shown  that  the  model  can 
successfully  predict  the  nonlinear  behavior  of  concrete  under  these 
complex  loading  conditions.  This  appears  to  provide  strong  indication 
that  the  present  model  has  captured  the  basic  constitutive 
characteristics  of  plain  concrete  quite  well. 


It  should  be  noted  that  the  scope  of  the  present  model  is 
currently  limited  to  stress  levels  in  plain  concrete  that  do  not 
produce  significant  cracking.  The  model  is  isotropic  and  therefore 
incapable  of  describing  stress- i nduced  anisotropy  produced  by 
significant  microcracking  with  preferred  orientation.  In  many 
practical  applications,  this  limitation  may  not  prove  to  be  overly 


SECTION  2 


THERMODYNAMIC  FOUNDATIONS  OF  THEORY 

Endochronic  plasticity  is  founded  on  the  theory  of  irreversible 
thermodynamics  of  internal  variables.  The  foundations  have  been  laid 
in  detail  elsewhere  (see  Ref.  21,  for  example)  and  will  not  elaborated 
upon  here.  For  the  sake  of  completeness,  however,  we  begin  by 
recalling  the  basic  equations. 

2.1  THERMODYNAMIC  RELATIONS. 

It  is  stipulated  at  the  outset  that  concrete  is  macroscop i ca I  I y 
homogeneous  and  isotropic  even  though  at  the  local  level  it  is  highly 
heterogeneous  and  anisotropic.  Its  mechanical  response  to  strain  (or 
stress)  is  also  highly  dissipative  and  therefore  inelastic.  In  the 
thermodynamic  formulation,  interior  dissipative  mechanisms  are 
formally  represented  by  internal  variables.  Isotropy  permits  the 
separation  of  such  mechanisms  into  deviatoric  and  hydrostatic  ones, 
represented  respectively  by  the  internal  variables  qr  and  pr,  where  qr 
are  tensors  of  second  order  and  pr  are  scalars.  In  the  context  of 
small  deformation  theory,  f,  the  free  energy  of  deformation,  is 
quadratic  in  the  strain  tensor  and  the  internal  variables  and  can 
therefore  be  partitioned  in  the  case  of  isotropic  materials  into  its 
deviatoric  and  hydrostatic  parts.  Thus 
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where  e  is  the  deviatoric  part  of  the  strain  tensor  £  and  e  is  the 
hydrostatic  part.  Thus 
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where  I  is  the  unit  matrix.  It  has  been  shown  that  in  the  case  of 
materials  which  behave  elastically  at  the  onset  of  unloading  or 
reloading  (points  of  reversal  in  the  sign  of  the  strain  rate)  can 
be  further  partitioned  into  an  elastic  and  plastic  part  such  that 


s  is  the  deviatoric  stress  tensor  and  u  is  the  elastic  shear  modulus. 
~ 

In  a  simi lar  fashion 


(2.9) 


where 


(2.10) 

(2.11) 


and 


(2.12) 

(2.13) 


a  is  the  hydrostatic  stress  and  K  is  the  elastic  bulk  modulus.  In 
para  I  lei  to  Eq.  (2.4) 

2  =  §  ♦  1°  ■  (2  14) 

The  topological  interpretation  of  the  above  partitioning  in  the  devia¬ 
toric  case  is  given  in  Figure  1.  The  same  i ntepretat i on  applies  to 
the  hydrostatic  case. 
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Figure  1.  An  elastic  system  and  a  plastic  system 
series . 


Figure  2.  Elastic  response  at  unloading  and  reloading. 


The  thermodynamic  relations  that  apply  are  the  fol lowing: ^  > 


ai;d9r  >°  .  Ild9rll  <° 


a  = 


*3 

3ef 


ar  dpr  >  0 
rr 


dpr|  >  0 


Also 


(2.15) 


(2.16) 


(217) 


(2.18) 


(2.19) 


(2 . 20) 


The  above  equations  are  obtained  by  applying  the  theory  of  irrever¬ 
sible  thermodynamics  to  each  individual  system  in  Figure  1  Strict 
interpretation  of  the  topological  arrangement  suggests  that  cannot 
depend  on  e,e  and  conversely  ^  cannot*  depend  on  ep. 


*.  However,  in  the  actual  material,  internal  changes  associated  with 
the  plastic  strain  may  affect  the  elastic  moduli  fie  and  «e  and  thus 
render  and  ^  dependent  of  gp  and  ep. 


We  note  again  that  the  free  energies  all  have  quadratic  forms 
which  have  been  shown  to  be  reducible  to  the  canonical  forms  given 
below 


*1 


2 


(2  21) 


"h  -  5  V£ 


f  el  2 


(2.22) 

(2.23) 

(2  24) 


It  is  physically  appealing  to  introduce  at  this  juncture  the 
concept  of  internal  forces  Qr  and  Pr  such  that 


(2.25) 


(2.26) 


It  has  been  shown  in  previous  work  that  qr  and  pr  represent  internal 
displacements  either  in  terms  of  assemblies  of  mechanical  mode  I  s^^, 
slip  in  crystals^*)  or  motion  of  atoms  grouped  together  in  terms  of 
specific  energy  states  .  Then  in  view  of  Eqs .  (2.16)  and  (2.25)  on 
one  hand  and  Eqs.  (2.17)  and  (2.26)  on  the  other, 
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These  may  also  be  given  in  terms  of  the  internal  forces  and  Pr  upon 
use  of  Eqs.  (2.25)  and  (2.26).  Thus  when  nr  and  Kr  are  constant, 


nr  dSr  ds 

^(1)  ^r  *  dt  “  ^r  dt 


(2.38) 


K  dP  .  o 

_L_  p  +  _! 1  _  K  dll 

^(2)  rr  dt  ~  r  dt 


(2  39) 


2.3  THE  NOTION  OF  INTRINSIC  TIME. 


In  materials  which  are  history  dependent  but  rate- i nd i f f erent 
"newton i an"  time  is  not  a  proper  time  scale  for  the  formulation  of 
evolution  equations.  The  difficulty  lies  in  t*-e  fact  that  in  such 
materials  the  stress  is  a  function  of  the  deformation  path  but  not  of 
the  rate  at  which  the  path  is  traversed.  This  fact  was  the  basis  of 
the  idea  proposed  by  the  first  author  in  1971^®)  that  every  material 


osses  an  intrinsic  time  scale  with  respect  to  which  the  evolution 


We  stipulate  that  the  tensor  P  is  isotropic, 
representati on 


Thus  it  admits  the 


P i j k£.=  kl5ijfik£  *  k2  5ikfij£. 


(2.41) 


which  when  combined  with  Eq.  (2.40)  admits  the  result 


2  -  >,2  aS 


6S  =  k l  de-k  d£P.  .  1.2  d£Pk  d£p_ 


(2.42) 


or,  within  an  immaterial  multiplicative  constant, 


~2  IUpII2  +  k2  |deP|2 


dz  =  de 


(2.43) 


dep  2  =  de?de?.,  I de ^ I  is  the  absolute  value  of  dep  and  k 

ij  i  j 


where 

is  a  material  parameter. 


If  one  retains  the  form  of  the  evolution  Eqs.  (2.31)  and  (2.32) 
and,  therefore,  (2.34)  and  (2.35)  and  regards  dz  in  Eq.  (2.40)  as  an 
intrinsic  time  measure  then  in  terms  of  this  time  scale  z  Eqs.  (2.34) 
and  (2.35)  become 


(2.44) 


(2.45) 
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These  two  Eqs.  lead  to  the  counterparts  of  Eqs .  (2. 38)  and  (2.39)  as 
shown  below. 


dA.  d± 

, (1)  *r  +  dz  ~  ^r  dz 


K  dP  .  p 

_ r_  p  ^  _ r  _ 

.  (2)  r  dz  ~  r  dz 

r 


(2.46) 


(2  47) 


Equations  (2.15)  and  (2.17)  give  the  relations  between  the  stresses 
and  the  internal  forces.  Specifically  in  view  of  Eqs.  (2.22)  and 
(2.24) 


!=I)r 


»  =  JZ  p 

r 


(2  48) 


(2  49) 


Equations  (2.43)  and  (2.46  -  2.49),  together  with  Eqs.  (2.8)  and 
(2.12),  provide  a  complete  constitutive  description  of  a  dissipative, 
rate- i ndependent  (plastic)  material,  subject  to  the  above 
st i pu I  at i ons . 

2.4  INTEGRAL  FORM  OF  THE  EQUATIONS 

In  general  the  resistance  tensors  bj;^  and  are  not  constant 
but  depend  on  other  variables  associated  with  the  deformed  state  of  a 
material.  A  propos  of  this  eventuality  we  set 


where  F,.  and  F^  are  specified  on  the  basis  of  the  observed  response  of 
a  particular  material.  This  question  will  be  dealt  with  in  the  next 
section.  Let 


where  b 


(1) 

or 


and  b 


(2) 

or 


are  constants. 


(2.52) 


(2.53) 


Further  we  define  intrinsic  time  scales  zs  and  z^  by  the  relation 


.  dz 

dz  =  c~ 
s  F 

(2. Ska 

s 

.  dz 

H  *  kFH 

(2.54b 

where  k  in  Eq.  (2.54)  is  the  same 
definition  of  intrinsic  time  given 
in  Eq.  (2.54)  is  not  purely  formal 
hydrostatic  response  so  derived  is 
later. 


as  the  parameter  k  appearing  in  the 
by  Eq.  (2.43).  The  appearance  of  k 
but  has  the  advantage  that  the  pure 
independent  of  k,  as  will  be  shown 


& 


In  the  light  of  the  above  definitions,  Eqs.  (2.46)  and  (2.47) 


become 


dSr  dep 

‘rlr'dT^rF 
s  s 


(2.55) 


dP 


X  P  ♦  ~  =  K  , 
r  r  dz^  r  dz^ 


(2.56) 


Let  the  reference  configuration  be  one  at  which  z^  =  zs  =  0.  We 
define  a  "natural  configuration"  as  one  at  which  JJr(0)  =  0,  Pr  =  0. 
There  are,  of  course,  configurations  (initially  strained  or  arising 
because  of  a  nonequilibrium  thermomechanical  reference  state)  in  which 
Qr(0)  /  0,  Pr(0)  t  0.  Thus  for  greater  generality  we  set 


9r(0)  =  9,  ,  Pr(0)  =  p° 


(2 . 57a , b) 


Equations  (2.55)  and  (2.56)  may  now  be  integrated  subjected  to  the 
initial  conditions  (2.57a,b)  to  give 


-a  z 

i  n°  r  s 

!r  ■  9r«  ♦  M, 


s  -«  (z  -2)  dep 

•  5F  d2 


(2.58a) 


-X  zu 

P  =  P°e  r  H  ♦  K 
r  r  r 


-Xr<Vl'>  dep 


dz 


r  dz' 


(2.58b) 
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*  * 

■M 


A 


S 


3 


si 


> 

V 

V 

>, 

> 

> 


-  •  ^ 


Use  of  Eqs.  (2.48)  and  (2.49)  gives  the  stress  in  terms  of  convolution 
integrals  of  the  plastic  strain  history  as  follows 


5  =  Z  9° 


-a  z 

o  r  s 
e  ♦ 


3ep 

'(V2')  ar-dz' 


(2.59) 


v—  o  VZH  ,N  ,  . 

^  =  Zpr  e  ♦  *(2H‘Z  )w  dz 

r  Jo 


(2 . 60) 


where 


p  =  E  r 


(2.61) 


*  =  lKre 


^rZH 


(2.62) 


The  condition  that  s  =  0  in  the  reference  state  requires  that 


Z  9?  =  2 


(2.63) 


whereas  the  condition  that  a  =  0  i n  the  refrence  state  requires  that 


X>°  =  0 


(2.64) 


If  the  reference  state  is  a  "natural"  state,  then 
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a  = 


'(V2')  dT  d2' 


(2  66) 


2.5  SINGULAR  KERNELS. 

It  was  shown  by  the  first  author^*)  that  if  the  kernel  p(zs)  is 
weakly  singular  in  the  sense  that  p( 0)  =  *  but  integrable, 

2s 

/)[zs]dzs  <  «  ,  zs<  «  (2.67) 

.  0 


then  the  hysteresis  loop  of  a  material  in  the  first  quadrant  of  the 
stress-strain  space  is  closed.  It  was  also  shown  that  the  material 
response  at  unloading  and  reloading  points  is  elastic  in  the  sense  of 
Figure  2.  The  reason  for  this  is  that 


dr 

1  -  dT 

dr 

d7P 

CL 

“O 

< 

B  d7p 

(2.68) 


because  of  the  singularity  of  the  kernel.  The  elastic  response  at 
these  points  is  then  the  response  of  the  elastic  system  of  Figure  1. 


In  numerical  computations  the  singular  kernel  is  approximated  by 
a  Prony  series  as  in  Eq.  (2.61),  in  which  care  is  taken  so  that  p(Q), 

where 


P( 0)  =  Z  Pr  ,  (269) 

r 

is  a  very  large  number  in  the  sense  that 

P(0)  »  pe  (2.70) 

where  is  the  elastic  slope  at  the  loading  point  A  of  Figure  2. 
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SECTION  3 


SPECIFIC  FORM  OF  THEORY  FOR  PLAIN  CONCRETE 

The  application  of  endochronic  plasticity  to  metals  is  facili¬ 
tated  by  three  assumptions  which  are  quite  realistic: 

(i)  Under  moderate  hydrostatic  stress,*  the  hydrostatic 
response  of  metals  is  elastic. 

(ii)  A  constant  moderate  hydrostatic  stress  does  not  affect 
the  mechanical  response  of  a  metal  in  shear. 

(iii)  Shearing  at  constant  hydrostatic  stress  does  not 
induce  a  change  in  the  hydrostatic  strain. 

In  concrete  the  above  assumptions  are  not  realistic. 

In  reference  to  item  (i)  above,  the  hydrostatic  behavior  of 
concrete  in  compression  is  illustrated  in  Figure  3.  It  is  initially 
convex  becoming  concave  and  asymptotically  elastic.  Upon  unloading  at 
A  to  a  point  B  a  significant  amount  of  plastic  strain,  eg  results. 
There  is  obviously  a  great  deal  of  hardening  taking  place,  which 
affects  dramatically  the  subsequent  I oad i ng-un I oad i ng-re I  cad i ng 
behavior  as  illustrated  in  Figure  3.  While  in  metals  hardening,  in 
general,  is  the  result  of  multiplication  of  dislocations,  in  concrete 
it  is  very  much  a  function  of  compaction.  In  both  cases  however  the 
agents  of  hardening  are  the  resistance  coefficients  br  and  in  this 
particular  instance,  .  Thus  while  in  metals  one  achieves 


•.  Of  the  order  of  the  yield  stress  in  tension. 
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hardening  by  letting  br  be  an  increasing  function  of  2,  in  concrete, 
to  account  for  the  compaction  effect  in  the  hydrostat,  we  have  set 


b(2)  =  b(2)  FH[fP] 


(3.1) 


where  is  a  monoton i ca I  I y  increasing  function  of  gP. 

In  reference  to  item  (ii),  the  shear  response  depends  strongly 
on  the  existing  level  of  hydrostatic  stress.  This  effect  is  illus¬ 
trated  schematically  in  Figure  4  where  the  shear  stress  is  plotted 
versus  the  shear  strain  under  monotonic  loading  conditions  and  at 
various  levels  of  constant  hydrostatic  stress.  In  regard  to  item 
(iii),  the  mechanical  response  of  concrete  shows  strong  shear-hydro¬ 
static  interaction  in  that  shearing  at  constant  hydrostatic  stress 
produces  a  significant  change  in  the  hydrostatic  strain  and  vice 
versa.  As  will  be  shown  subsequently,  endochronic  plasticity  accounts 
for  this  effect  through  the  intrinsic  time  and  specifically  by  virtue 
of  the  coupling  constant  k  which  appears  in  Eq.  (2.43).  Note  that  in 
materials  which  are  plastically  incompressible,  de^  =  0.  In  such 
materials  shear-hydrostatic  interaction  is  absent. 

3.1  HYDROSTATIC  RESPONSE. 

Determination  of  the  material  functions  ^(z^)  and  F^(eP) .  It  is 
assumed  at  the  outset  that  the  initial  state  is  a  natural  state  in 
which  case  the  hydrostatic  response  is  given  by  Eq.  (2.66)  which  we 
write  in  the  explicit  form 


■*>>>'* 


Under  monotonic  loading  conditions,  and  in  view  of  Eqs.  (2.43)  and 
(2.54),  it  follows  that 


d£ 

dzH 


FH 


Thus  letting 

d?H  =  ldepl 


it  fol lows  that 


k  dep  _  d?H 
H  kFH(ep)  FH(epl 


dzu  = 


Under  monotonic  straining  conditions  dzH  =  deP  hence 


d?H  = 


d6f 


FH(ep) 


Thus  Eq.  (3.2)  becomes 


z')FH(ep')dz' 


Two  material  functions,  therefore,  are  necessary 
determine  the  hydrostatic  response:  the  kernel 

hardening  function  F^(ep) . 


(3.3) 

(3.4) 

(35) 

(36) 

(3.7) 

and  sufficient  to 
^(z^j)  and  the 


3.1  1  The  Hardening  Function  F^. 


This  function  will  be  examined  in  reference  to  the  hydrostatic 
stress-strain  curve  shown  in  Figure  5.  The  curve  ON  represents  the 
stress  response  corresponding  to  F|^  =  1 .  In  the  initial  stage 
(eP  <  OA)  F^  is  substantially  linear.  In  the  stage  OA  <  ep  i  OB  it  is 
substantially  hyperbolic.  In  the  final  stage  eP  >  OB  a  ■*  *  as 
eP  -►  eP(  since  there  is  a  limiting  material  compressibility  beyond 
which  the  material  cannot  be  compressed  no  matter  how  high  a  stress  a 
i s  app lied. 

We  shall  take  the  position  that 

Fh  =  (3.8) 

which,  in  the  range  cp  <  OA,  we  can  write  as 

FH  a  1  ♦  J£P  ♦  o(eP)2  (3.9) 


The  linear  representation  of  F^  inherent  in  Eq.  (3.9)  affords  an 
explicit  method  of  determination  of  the  value  of  P  and  the  kernel 

^(zh) ■ 

3.1.2  Method  of  Analysis. 


(3.10) 


It  follows  from  Eqs .  (3.67)  and  (3.9)  that 


/yy  /«* 


s 


fz 

s 

P[Zs 

.  0 


(3 . 20) 


where 


dz  = 


de^ 


!(deP): 


(3.21a) 


and 


dzs  =  p  (3.21b) 

s 

in  accordance  with  Eqs.  (2.42)  and  (2.54).  Thus  two  material  func¬ 
tions,  />(zs)  and  Fs  and  the  constant  k  are  sufficient  to  determine  the 
shear  stress  from  the  plastic  strain  history.  It  has  been  found 
experimentally  that  the  failure  stress  of  concrete  in  the  T-plane 
depends  on  the  prevailing  hydrostatic  stress  and  the  direction  of 
loading  in  the  f-plane.  To  account  for  these  two  effects  it  is  neces¬ 
sary  that  Fs  depend  on  a  and  the  Lode  angle  6.  Thus  we  set 

Fs  =  Fs(<7,3)  (3.22) 

subject  to  the  normalization  condition  Fs(0,0)  =  1. 

3.2.1  The  Kernel  p(zs) . 

The  form  of  the  function  p( zs)  is  found  most  conveniently  by 
means  of  a  pure  shear  test  in  the  absence  of  hydrostatic  stress.  In 
this  case 
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If  the  non-vanishing  stress  component  is  t  and  the  corresponding 
plastic  strain  component  is  7P,  then  in  view  of  Eq.  (3.21)  we  can 


wr  i  te 


dz  =  f2!d7pl 


dz  =  dz 
s 


Thus  using  Eq.  (3.20) 


(3.23) 


T  =  —  m(J27P) 

(2 


(3.24) 


where 


M(x)  =  p(z')dz 


for  all  f i ni te  x .  Thus 


(3  25) 


=  a-  =  p(n7P) 

<f(j27p)  d7p 


(3.26) 


Thus,  if  one  plots  P2t  versus  l27P,  the  slope  of  the  curve  is  the 
kernel  function  a(x)  where  x  =  1"27P. 


V 


3.2.2  The  Function  Fg. 


The  determination  of  the  form  of  Fg  and  the  value  of  k  is  more 
complex  due  to  the  strong  coupling  between  the  hydrostatic  and  dev  i  - 
atoric  responses  brought  about  by  the  presence  of  k  in  Eq.  (3.21). 

Consider  for  instance  a  shear  stress  field  in  the  presence  of  a 
constant  hydrostatic  stress.  More  precisely  a  hydrostatic  stress 
field  is  applied  until  a  -  0^,  whereupon  a  is  held  constant  at  ffj  and 
a  shear  stress  field  f  is  then  applied,  so  that  while  r  is  being 
applied  the  total  stress  field  is  of  the  form 


l 


where  H(z)  is  the  Heaviside  step  function.  Let  z^  =  at  the  ter¬ 
mination  of  the  hydrostatic  test,  i.e.,  at  a  =  ffj.  Then  in  view  of 
Eq.  (3.28) 


o 

H 


ZH)^dzH 


since 


da 

dz. 


=  0  when  z^  >  z 


H 


(3.30) 


Since  the  application  of  shear  stress  will  cause  a  change  in  the 
shear  strain,  z  will  change  by  virtue  of  Eq.  (3.21),  and  so  will  z^  by 
virtue  of  Eq.  (2.54).  Thus  the  integral  on  the  righthand  side  of  Eq. 
(3.30)  will  change  and  hence  so  will  eP.  Thus,  in  the  presence  of 
hydrostatic  stress,  dz  is  no  longer  given  by  Eq.  (3.23)  but  by  Eq. 
(3.21),  which  is  influenced  by  the  value  of  k. 

The  problem  simplifies  somewhat  and  becomes  amenable  to  analysis 
if  one  adopts  a  simpler  form  of  the  hydrostat.  With  reference  to  Eq. 
(3.2),  the  kernel  0(z^)  decays  with  sufficient  rapidity  to  be  approxi¬ 
mated  by  a  Direc  5  function.  Thus  we  set 


*[zh)  ■  V(zhJ 


in  which  case  by  virtue  of  Eq .  (3.2),  it  follows  that 


a 


=  K 


de 
0  dz 


P 

H 


(3.31) 


(3.32) 


Vvl 


Specifically  under  purely  hydrostatic  monotonic  loading  conditions  Eq. 
(3.3)  applies  in  which  case 


a  =  K  (l  ♦  ptp) 
o  r 


(3.33) 


by  virtue  of  Eqs.  (3.9)  and  (2.43).  Equation  (3.33)  is  then  basically 
that  of  a  plastic  linearly  hardening  model.  In  particular  if  aj  is  on 
the  linear  part  of  the  hydrostatic  stress-strain  curve,  Eqs.  (3.33)  is 
certainly  adequate. 

It  follows  from  Eqs.  (3.32)  that 


dfP  =  Kq  dZH 


(3.34) 


where  ffj  is  the  (constant)  value  of  the  hydrostatic  stress  during 
shear i ng. 

Adopting  the  sign  convention  commonly  used  in  civil  engineering 
whereby  a  and  are  positive  when  compressive,  we  note,  in  view  of 
Eq.  (3.34),  that  de^  will  be  positive,  and  hence  the  compressive 
strain  will  increase,  during  shearing. 


In  view  of  the  above  observations  we  introduce  the  following 


notation 


dc  =  deP  ,  dcu  =  ldep 


(3 . 35a  .  b) 


I. 

V 


Thus,  because  of  Eq.  (2.43), 


,  2  .2  .2,2 

dz  =  d£s  ♦  k  d^H 


(3.36) 


We  now  combine  Eqs .  (2.54)  and  (3.32)  to  obtain 


«  -  kKo(]  •  ^zh)  ST 


(3  37) 


where  use  was  made  of  Eq.  (3.33),  (3.35b)  and  the  fact  that  during  the 
hydrostatic  as  well  as  the  shearing  processes  dep  >  0.  Equations 
(3.36)  and  (3.37)  combine  to  give  the  result 


dz 


1 

1*1 

1  - 

Ko(1  ♦  *h) 

1/2 


* 


(3.38) 


Now  let  fH(=e^)  *3e  value  of  €p  at  the  completion  of  the  purely 

hydrostatic  loading  so  that 


(3 . 39) 


where  is  the  hydrostatic  strain  induced  by  shear  loading  at  cons¬ 
tant  pressure.  Also  let 


where  z'  is  the  value  of  z  at  the  completion  of  the  purely  hydrostatic 
loading.  Without  going  into  extensive  analytical  detail  Eqs.  (3  36), 
(3.37)  and  (3.38)  combine  to  give  the  following  results: 


(3.41) 


where 


(3  42) 


(3.43) 


where 


x 


1  ♦ 


V 


(3  44) 


Relations  (3.41)  and  (3.43)  will  prove  useful  in  determining  the  form 
of  Fs  and  the  value  of  k. 

We  now  return  to  the  question  of  the  shear  stress  response  to 
increasing  shear  strain  in  the  presence  of  constant  pressure,  i.e.,  to 
the  case  where  the  stress  field  is  given  by  Eq.  (3.27).  In  view  of 
Eq.  (3.20),  we  can  write 


However,  7P  =  0  during  pure  hydrostatic  loading,  i.e.,  i 
0  <  zs  <  zj.  Thus,  if  we  set  w  =  zs  -  z^  and  recall  Eq. 
fol lows  that 


r 


»  d<vp 

P(w  “  d" 

0 


We  note  that 


d7P  d7P  d?s  _  1  d?s 

dw'  ~  dfs  dw  dw 


in  view  of  Eq.  (3.35a).  Therefore,  as  a  result  of  Eqs. 
(3.41) 


F 

s 


2aFsw 
1  ♦  2arsw 


1/2 


»here  in  this  case 


ps  =  Fs(o,0) 


i.e.,  F^  is  constant  during  the  shearing  process. 


the  range 
(3.45),  it 

(3.45a) 

(3.46) 

3.21b)  and 

(3.47) 

(3  48) 


Thus  as  a  result  of  Eqs.  (3.45a),  (3.46)  and  (3.47) 


*:»1 
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rv-  -  'U  ♦  2aFwl 
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It  is  shown  in  Appendix  I  that 

i ■  fW  2aF  z ' 

I im  -  s 

z-*«  ^'z  z  '  1  *  2aF  2 

.0  S 


(3 . 50) 


where  M(z)  is  given  by  Eq.  (3.25)  and  Mw  =  M(»)  <  ®.  Hence,  denoting 
the  asymptotic  value  of  the  shear  stress  by  r,,  it  follows  from  Eqs. 
(3.49)  and  (3.50)  that 


T  =  —  M 
"  /2  * 


During  shearing  at  zero  pressure  one  may  set  F  =1.  Thi 


(3.51) 


T®(glJ  - 

r«(0)  * 


(3.52) 


Equation  (3.52)  determines  the  form  of  Fs  under  constant  9.  Keeping  a 
constant  and  performing  shearing  tests  each  at  a  different  value  of  8 
one  is  able  to  determine  the  form  of  functional  dependence  of  Fg  on  8. 

3.2.3  Determination  of  the  Kernel  ^(w)  in  the  Presence  of  Pressure. 

In  the  course  of  the  development  of  the  theory  in  conjunction 
with  the  experimental  data  of  Scavuzzo,  et  aL.^^)  ; t  transpired  that 


all  shear  experiments  were  performed  in  the  presence  of  pressure  i.e., 
under  conditions  where  t  0.  In  this  event  Eq.  (3.24)  no  longer 
applies  and  Eq.  (3.49)  must  be  used  for  the  determination  of  p(w) .  To 
this  end  we  set  Fs  =  1  at  a  reference  pressure  ffp,  which  is  the  one  we 
use  in  the  determination  of  /j(w),  and  which  must  lie  on  the  straight 
line  segment  of  the  hydrostat  to  ensure  that  Eq.  (3.49)  is  applicable. 
Thus  setting  Fs  =  1  in  Eq.  (3.49)  we  obtain  the  following  relation  for 


T  =  7=  ply  -  y')(t Tur)  '  d»' 
w  Jo 


(3.53) 


where  as  before 


(3.54) 


and  since  y  =  w  when  Fs  =  1 . 

Thus  given  r  (which  is  measured  experimentally)  and  a,  the  task 
at  hand  is  to  solve  the  integral  equation  (3.53)  and  thus  obtain  the 
functional  form  of  p. 

Before  we  proceed  with  the  solution  we  take  care  to  note  that  r 
is  measured  as  a  function  of  the  shear  strain  7  and  is  therefore 
reducible  to  a  function  of  the  plastic  strain  7P .  In  the  presence  of 
monotonic  shearing  conditions  and  in  view  of  Eq.  (3.35a) 


dC  =  J2  d7P 


(3.55) 


Furthermore,  since  at  the  initiation  of  the  shear  test  gs  =  0,  it 
fol I ows  that 


(3.56) 


Thus  in  view  of  the  relation  between  £s  and  y  (Eq.  3.44),  £s  can 
be  expressed  as  a  function  of  y.  Thus  the  function  on  the  lefthand 
side  of  Eq.  (3.53)  is  a  known  function  of  y. 


3.2.4  Numerical  Procedure  for  Determining  the  Kernel  Function  />(zs). 
We  proceed  to  give  a  numerical  solution  to  the  equation 


r(y) 


y 

p(y  -  y')  g(y')dy' 
o 


(3.57) 


where  the  object  is  to  determine  p(y)  given  the  functions  r(y)  and 
g(y) .  Note  that  Eq.  (3.53)  is  a  particular  case  of  Eq.  (3.57)  where 


9(y)  =  f2  (r? 


2ajL_) 


1/2 


2ayJ 


(3.58) 


Let  yn  be  a  generic  value  of  the  variable  y.  We  divide  the  interval 
[0,y]  into  suitably  small  equal  subintervals  &  such  that 


The  integral  on  the  righthand  side  of  Eq.  (3.57)  can  then  be  written 
in  the  form 


r 


r 

n 


'A 

p[yn  -  y']g(y')dy'  ♦  •  •  • 
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'nA 

p[ yn  -  y')g(y')dy' 

(n-l)A 


or 


rn  =  21  gr{w[(n  -  r  *  1) A]  -  M[(n  -  r)A]} 
r=l 

where  y  virtue  of  the  mean  value  theorem 

gr  =  g[y?j  . 

(r  -  1)A  <  y»  <  rA 


and 


M(rA) 


'rA 

P(  y')dy 

o 


The  quantity  gr  may  be  determined  to  the  required  degree  of 
by  making  the  interval  A  sufficiently  small. 

In  consequence,  we  have  the  relations 


(3.60) 


(3.61) 

(3.62) 

(3.63) 

(3.64) 

accuracy 
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(3.65) 


r2  =  (2A)  ♦  Ag2M(A) 


rn  =  AgjM(nA)  ♦  Ag2M[(n  -  1) A] 


♦  ...  +  Ag  M (A) 


where 


Agr  =  gr  -  gr_: 


(3.66) 


This  is  a  system  of  a  linear  simultaneous  equations  in  M(A),  M(2A)  ... 
M(nA) .  Since  the  Agr  are  known,  then  knowing  Tj,  t2  ...  rn  allows  one 
to  determine  M(rA),  r  =  1,2  ...  n,  by  successive  substitution.  The 
following  a  I gor i thm  app I ies : 


If 


T 

n 


m=l 


(3.67) 


then 


M(rA) 


(3.68) 


Knowledge  of  M(y)  thus  allows  the  determination  of  p( y)  since 


My)  = 


(3.69) 
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SECTION  4 

EVALUATION  OF  MODEL  PARAMETERS  FROM  PLAIN  CONCRETE  DATA 


In  the  preceding  section,  a  specific  form  of  the  general  endo- 
chronic  model  described  in  Section  2  was  given  for  plain  concrete,  and 
a  basic  approach  was  given  for  determining  the  material  constants  in 
the  model  from  experimental  data.  In  this  section,  the  model  is 
applied  specifically  to  the  data  on  plain  concrete  recently  reported 
by  Scavuzzo,  et  aL.(^)  The  model  is  first  fit  to  a  small  subset  of 
the  data,  and  the  resulting  model  is  then  proof-tested  against  the 
data  from  a  wide  variety  of  the  complex  strain  path  tests  not  used  in 
fitting  the  model.  The  predictive  capabilities  of  the  model  are 
demonstrated . 

4.1  DETERMINATION  OF  MODEL  PARAMETERS. 

Only  a  small  subset  of  the  data  reported  by  Scavuzzo,  et 
is  needed  to  determine  the  values  of  the  model  parameters;  for  this 
purpose,  pure  hydrostatic  compression  data,  data  on  the  response  to 
shear  at  several  different  fixed  hydrostatic  pressures  and  triaxial 
failure  data  at  various  confining  pressures  proves  sufficient. 
Details  of  the  manner  in  which  this  can  be  accomplished  are  given 
below. 

4.1.1  Hydrostatic  Parameters. 

The  most  extensive  virgin  pure  hydrostatic  compression  curve 
reported  by  Scavuzzo,  et  a is  depicted  by  the  dotted  line  in 
Figure  7,  covering  the  range  of  pressures  up  to  nearly  14  ksi.  Shown 
also  in  Figure  7  by  the  solid  line  is  the  virgin  hydrostatic 
compression  curve  actually  used  in  determining  the  model  parameters. 


An  examination  of  the  virgin  pure  hydrostatic  compression  legs  of  many 
of  the  complex  strain  path  tests  revealed  that  the  dotted  curve  in 
Figure  7  was  probably  too  soft.  Therefore,  from  this  examination,  the 
curve  shown  by  the  solid  line  was  constructed  and  appears  to  provide  a 
more  accurate  portrayal  of  the  normal  virgin  hydrostatic  behavior;  the 
linearization  of  the  upper  portion  of  the  hydrostatic  curve  was  felt 
to  be  justified  in  view  of  the  fact  that  none  of  the  complex  stress 
path  tests  reported  by  Scavuzzo,  et  aL.  exceeded  hydrostatic  pressures 
of  12  ksi .  Note  that  the  curves  shown  in  Figure  7  are  expressed  in 
terms  of  octahedral  normal  stress  aQ  versus  octahedral  normal  strain 
eQ.  Since  aQ  =  1/3^7 j  j  and  eQ  =  l/3Cjj,  the  mean  stress  (pressure)  a 
is  synonymous  with  ffQ,  and  the  volumetric  strain  £  is  equal  to  3  £0. 

Based  on  the  curves  given  in  Figure  7,  and  in  agreement  with 
Stankowski  and  Gerstle,^®)  the  bulk  modulus  K  was  assigned  the  value 


K  =  2.1  x  10°  ksi 


(4.1) 


Using  the  solid  curve  in  Figure  7,  together  with  the  expression 


ep  =  e  -  - 
e  -  e  K 


(4.2) 


and  the  value  of  K  given  by  Eq.  (4.10),  a  curve  of  a  versus  was 
obtained  which  is  shown  in  Figure  8.  From  the  straight  line  portion 
of  this  figure,  one  finds  by  recalling  Eq.  (3.18)  given  earlier  that 


Using  these  values  of  p  and  aQ,  a  plot  was  made  of  aQ  -  a*  versus  z^, 
using  the  curve  from  Figure  8,  and  the  definitions  of  C7«  and  z^  given 
earlier  by  Eqs.  (3.11)  and  (3.16),  respectively.  Six  points,  at 
equally  spaced  values  of  z^,  were  chosen  from  this  plot,  and  a  two- 
term  series  of  decaying  exponentials  of  the  form 


a  -  a * 
o 


.  t  C 

rtl  r 


was  fit  to  these  points  by  Prony’s  method.  This  led  to  the  following 
values  of  Cr  and  f)r: 


Cj  =  2.434  ksi 


Px  =  635 


C2  =  2.566  ksi 


fi2  =  2,289 


A  plot  of  Eq.  (4.4),  based  on  the  constants  given  by  Eq.  (4.5),  is 
shown  in  Figure  9,  where  the  six  points  to  which  Eq.  (4.4)  was  fit  are 


also  shown. 


By  comparing  Eqs.  (4.4)  and  (3.19),  it  then  follows  that 


C  = 

r  p  ♦  Xr 


Pr  =  P*  \ 


from  which  we  find 
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Figure  9.  Representation  of  the  oQ-a»  versus  z^  rela¬ 
tionship  by  a  two- term  series  of  decaying 

exponentials.  The  solid  circles  are  the 
equally-spaced  points  to  which  the  series  was 
fit. 


Figure  10.  Octahedral  normal  stress  versus  octahedral 
normal  strain  for  pure  hydrostatic  compression, 
including  unloading  and  reloading. 


(4.7) 


Kj  =  1.55  x  103  ksi,  \l  =  570 


K2  =  5.87  x  10°  ksi,  X2  =  2,224  , 


since  p  is  known  for  Eq.  (4.3).  Therefore,  Eq.  (3.14),  with  the  above 
values  of  the  constants,  completely  defines  the  hydrostatic  behavior 
of  the  model.  An  example  of  the  hydrostatic  behavior  predicted  by 
this  model,  including  unloading  and  reloading  response,  is  given  in 
Figure  10.  Comparison  of  this  figure  with  the  solid  curve  in  Figure  7 
shows  that  the  model  portrays  the  hydrostatic  behavior  quite  well. 

4.1.2  Shear  Parameters. 

As  noted  earlier  in  Section  3.2,  the  most  convenient  and  direct 
way  of  determining  the  shear  parameters  is  from  data  on  shear  response 
under  zero  hydrostatic  pressure.  Unfortunately,  Scavuzzo,  et  al .W) 
did  not  investigate  this  case  in  their  extensive  study.  They  did, 
however,  investigate  shear  behavior  under  several  different  fixed 
hydrostatic  pressures  in  Tests  4-9  to  4-11;  the  resulting  data  will  be 
used  to  determine  the  shear  parameters  and  the  shear- vo I umetr i c 
coupling  parameters  in  this  section. 

A  composite  of  the  results  for  those  portions  of  Tests  4-9  to  4- 
11  which  consisted  of  shear  at  constant  hydrostatic  pressure  is  given 
in  Figure  11.  An  inspection  of  this  figure  reveals  some  apparent 
discrepancies  in  the  data.  From  the  general  characteristics  of  plain 
concrete  reported  in  the  literature,  the  initial  slope  at  the  onset  of 
shear  generally  increases  with  confining  pressure,  and  the  shear 
stress-shear  strain  curves  can  be  expected  to  fall  above  each  other  — 
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not  intersect  —  as  the  confining  pressure  increases.  Tests  4-9  and 
4-10  follow  the  expected  pattern,  but  Test  4-11,  despite  being  at  the 
highest  confining  pressure  (8  ks i ) ,  does  not;  it  has  the  smallest 
initial  slope  and  exhibits  softer  response  than  Test  4-10.  There  is 
undoubtedly  scatter  in  these  data  but  since  the  tests  were  not 
repeated,  the  magnitude  of  the  data  scatter  is  unknown. 

In  Section  3.2,  an  approach  for  determining  the  shear  parameters 
from  shear  data  at  constant  hydrostatic  pressures  was  described,  which 
required  the  use  of  shear  data  obtained  at  fixed  hydrostatic  pressures 
that  lay  on  the  linear  portion  of  the  hydrostat.  While  the  test 
conditions  for  Test  4-11  meet  this  requirement,  i.e.,  shear  at 
aQ  =  8  ksi,  the  data  are  somewhat  suspect,  as  noted  above.  In  view  of 
this,  we  used  our  best  judgement  to  construct  a  curve  which  we  felt 
provided  a  more  realistic  description  of  the  response  of  plain 
concrete  to  shear  at  aQ  =  8  ksi,  and  this  curve  is  shown  by  the  solid 
line  in  Figure  11.  The  reference  pressure  ffp  was  therefore  set  at 
8  ksi,  and  the  shear  parameters  determined  from  the  solid  curve  in 
Figure  11 . 

The  shear  modulus,  2 fi,  and  its  variation  with  deformation  was 
also  determined  directly  from  the  curves  shown  in  Figure  11.  The 
shear  moduli  were  obta  ned  from  the  initial  slopes  of  the  loading 
curves,  as  well  as  from  the  slopes  of  the  initial  portions  of  the 
unloading  curves.  Attempts  were  made  to  correlate  the  resulting 
values  of  2p  against  the  plastic  volumetric  strain,  the  total 
volumetric  strain,  the  octahedral  shear  strain  and  the  plastic 
octahedral  shear  strain.  The  total  volumetric  strain  and  the  plastic 
volmetric  strain  provided  poor  correlation,  while  the  octahedral  shear 
strain  and  the  plastic  octahedral  shear  strain  led  to  reasonably  good 
correlations.  Adopting  the  plastic  octahedral  shear  strain,  a  linear 
fit  to  the  data  points  in  Figure  12  led  to  the  following  expression 


2)1  =  2Mo  -  m  7^ 


(4  8) 


in  which  2/i0  =  1.83  x  10®  ksi  and  m  =  1.42  x  10®  ksi  . 

By  using  Eq.  (4.8),  the  following  incremental  expression  for  the 
plastic  octahedral  shear  strain  7^ 


<  -  • 


(4.9) 


and  the  solid  curve  in  Figure  11,  the  dependence  of  rQ  on  7P  at  the 
reference  pressure  <7p  =  8  ksi  was  determined,  and  the  result  is 
depicted  in  Figure  13. 

At  this  point,  it  becomes  necessary  to  determine  the  shear-volu- 
metric  coupling  coefficient  k  and,  for  this  purpose,  we  consider 
Figure  14,  which  shows  the  manner  in  which  the  octahedral  normal 
strain  eQ  varies  with  the  octahedral  shear  stress  Tq  during  shear  at 
several  different  constant  hydrostatic  pressures.  In  the  absence  of 
contrary  information,  it  will  be  assumed  that  volumetric  changes  which 
take  place  during  shear  at  constant  pressure  are  completely  irrever¬ 
sible.  On  this  basis,  and  using  the  results  for  aQ  =  8  ksi  depicted 
in  Figures  13  and  14,  a  plot  of  the  plastic  volumetric  strain  ep 
versus  the  plastic  octahedral  shear  strain  7P  was  generated,  and  this 
is  shown  in  Figure  15. 


The  model  provides  a  theoretical  expression  for  the  dependence 
of  ep  on  7P  through  Eq.  (3.43),  given  earlier,  if  we  note  that 


( ksi) 


Figure  13.  Dependence  of  the  octahedral  shear  stress  rQ  on  the 
plastic  octahedral  shear  strain  for  a  fixed 
hydrostatic  pressure  of  8  ksi . 


Figure  14.  Dependence  of  the  octahedral  normal  strain  on  the 
octahedral  shear  stress  rQ  for  several  different 
values  of  constant  hydrostatic  pressure.  Data  are 
from  Tests  4-9  t  4-11. 


(4.10) 


ep  =  (x  -  1)  /  (a k) 

With  these  identifications,  the  dependence  of  on  7P  is  obtained  for 
a  specified  value  of  k,  since  the  constant  a  is  known  in  terms  of  k 
through  Eq.  (3.42).  In  this  manner,  plots  of  versus  7P  were 
generated  with  the  aid  of  Eq.  (3.43)  for  a  range  of  assumed  values  for 
k,  and  the  results  are  shown  in  Figure  15.  From  an  inspection  of  this 
figure,  the  curve  corresponding  to  k  s  1.5  appears  to  provide  the  best 
overall  description  of  the  data;  thus,  we  set 

k  =  1.5  (4.11) 


i  n  the  mode  I . 


In  Appendix  C,  the  following  expression  is  derived,  relating 
and  zs  during  shear  at  constant  pressure,  for  the  case  in  which  a  is 
at  the  reference  pressure  ffp  and  on  the  linear  portion  of  the  virgin 
hydrostat: 


Ss  =  k  (i2a  2s  I1  7  2a  zs  ’  ' °9 I | 2a  zs  -  l1  *  2a  ZS  I)  (412) 


During  shear  at  constant  pressure,  it  can  be  shown  that 


(4.13) 
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which,  when  combined  with  Eq.  (4.12),  provides  a  relationship  between 
7P  and  zs  during  shear  at  constant  pressure  for  a  prescribed  value  of 
a.  Adopting  the  value  of  k  given  in  Eq.  (4.11),  the  value  of  the 
constant  a  is  then  known,  and  the  resulting  expression  given  by  Eqs. 
(4.12)  and  (4.13)  between  7P  and  zs  is  as  shown  in  Figure  16.  Using 
this  result,  and  the  relationship  between  rQ  and  7P  given  in  Figure 
13,  the  dependence  of  rQ  on  zs  during  shear  at  constant  pressure  can 
be  constructed.  The  numerical  method  described  in  Section  3.2  may  now 
be  utilized  in  conjunction  with  the  rQ  versus  zs  relationship  just 
found  to  give  M(zs),  where  M  is  defined  by  Eq.  (3.64).  In  this 
manner,  the  result  shown  in  Figure  17  was  obtained.  Using  this  in 
conjunction  with  Eq.  (3.69),  the  shear  kernel  function  p(zs)  can  be 
found  by  differentiation  of  the  M(zs)  relationship.  A  two-term  series 
of  decay  exponentials  of  the  form 


'M  *  £  V 


-a  z 
r  s 


(4.14) 


was  fit  to  the  shear  kernel  by  Prony’s  method;  this  led  to  the  fol¬ 
lowing  values  for  the  Af  and  ar:* 


Aj  =  1.46  x  10°  ksi  ,  ffj  =  100 


A_  =  19.0  x  10°  ksi  ,  a  =  6,554. 


(4.15) 


*.  Because  the  numerical  inversion  of  Eq .  (3.53)  was  only  approximate 
(see  Figure  17),  the  p{z)  so  obtained  was  also  approximate.  As  a 
result,  an  adjustment  in  the  value  of  a,  to  that  listed  in  Eq. 
(4.15)  was  made  to  provide  a  better  description  of  the  shear 
behav i or . 


Finally,  on  the  basis  of  the  failure  data  given  by  Stankow- 
the  following  expression  for  the  function  Fg  was  adopted  for 
use  in  the  present  study: 


t*  ♦  B  a 

F  =  - ~~ 

s  t*  ♦  fi  a 


s  R 


(4.16a) 


in  which 


r»  =  2  ksi 


Ps  =  0.5  (4.16b) 

ap  =  8  ksi 

This  completes  the  specification  of  the  endochronic  concrete  model. 

4.2  INTERNAL  CONSISTENCY  OF  MODEL,  DATA  AND  FITTING  PROCEDURE. 

We  return  now  to  determine  how  well  the  plain  concrete  model 
just  specified  can  describe  the  data  to  which  it  was  fit,  namely,  the 
data  from  Tests  4-9  to  4-11  of  Reference  17.*  Such  an  exercise  does 
not,  of  course,  constitute  a  proof-test  of  the  model,  but  simply 
provides  a  test  of  the  internal  consistency  of  the  fitting  procedure. 
It  should  be  remembered  that  plain  concrete  data  exhibits  significant 
data  scatter  and,  in  comparing  a  model  prediction  against  an  indivi- 


A  comparison  between  the  measured  and  predicted  hydrostatic 
behavior  was  given  earlier  in  Figure  15,  and  will  not  be  repeated 
here. 
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dual  test  result,  one  cannot  expect  to  obtain  the  same  level  of  agree¬ 
ment  as  is  found  for  example  in  metals.  With  plain  concrete,  one  is 
looking  more  for  a  model  that  captures  the  overall,  general  character¬ 
istics  of  response. 


To  explore  the  ability  of  the  model  to  describe  the  data  from 
Tests  4-9  to  4-11,  the  model  was  driven  around  the  appropriate 
prescribed  stress  paths  for  these  tests,  using  a  stress-driven 
computer  program  based  upon  the  explicit  incremental  numerical  scheme 
described  in  Appendix  B.  The  results  from  these  calculations  are 
depicted  in  Figures  18  to  20.  As  an  inspection  of  these  figures 
reveals,  the  overall  agreement  between  predictions  and  the  data  is 
quite  good,  considering  the  uncertainty  in  the  data  discussed  earlier. 
The  figures  reveal  that  the  main  discrepancies  lie  in  the  unloading 
behavior,  with  the  model  showing  noticeably  more  curvature  of  the 
unloading  path  —  and  hence  more  hysteresis  —  than  the  data. 
Although  no  effort  was  made  in  the  present  study  to  refine  the  model 
so  that  it  would  more  closely  describe  the  unloading  behavior,  such  an 
improvement  is  straightforward  and  will  be  accomplished  in  a  future 


(a)  Octahedral  shear  stress  versus  octahedral  shear  strain 


Figure  19.  Measured  and  predicted  responses  to  shear  at  a 
fixed  hydrostatic  pressure  of  6  ksi .  (Test  4- 


(a)  Octahedral  shear  stress  versus  octahedral  shear  strain 


Figure  20.  Measured  and  predicted  responses  to  shear  at  a 
fixed  hydrostatic  pressure  of  8  ksi  (Test  4- 


SECTION  S 


PROOF-TESTS  OF  MODEL  OVER  COMPLEX  STRESS  PATHS 

In  the  preceding  section,  the  endochronic  constitutive  model 
developed  in  Section  3  was  applied  to  the  plain  concrete  data 
presented  in  Reference  17,  with  the  model  parameters  evaluated  from  a 
small  subset  of  the  reported  data.  The  ability  of  the  resulting  model 
to  accurately  portray  the  data  to  which  it  was  fit  was  also  described. 
The  test  program  reported  in  Reference  17  consisted  of  six  different 
series  of  stress-controlled  tests,  each  series  being  designed  to 
explore  a  particular  facet  of  plain  concrete  behavior.  In  that  which 
follows,  the  model  is  proof-tested  by  examining  its  capability  to 
predict  the  measured  response  around  a  selected  number  of  prescribed 
complex  stress  paths  from  this  test  series.  None  of  the  stress  path 
tests  considered  in  this  section  were  used  in  fitting  the  model 
parameters,  and  no  optimization  techniques  were  employed  to  achieve 
the  results  presented  below. 

5.1  SERIES  NO.  1. 

This  test  series  was  designed  to  explore  the  response  of  plain 
concrete  to  triaxial  load  cycles  which  do  not  exhibit  stress 
reversals.  Each  test  was  repeated  once.  The  model  was  applied  in  the 
present  study  to  Test  1-1  from  this  series,  for  which  the  prescribed 
load  paths  are  shown  in  Figure  21.  The  loading  consisted  first  of 
cyclic  hydrostatic  loading  up  to  8  ksi,  as  shown  in  Figure  21(a), 
followed  by  cyclic  deviatoric  loading  along  the  triaxial  compression 
path,  as  illustrated  in  Figure  21(b).  The  predicted  and  measured 
response  for  this  loading  history  are  given  in  Figure  22.  As  this 
figure  reveals,  the  predicted  and  measured  responses  are  in  good 
agreement,  considering  the  data  scatter  which  can  be  assessed  by 


(a)  Cyclic  hydrostatic  loading. 


(b)  Cyclic  deviatoric  loading  along 
the  tri axial  compression  path. 


Figure  21.  Stress  paths  for  Test  1-1. 


PRINCIPAL  STRAINS 


comparing  the  results  from  Test  1-1  with  those  from  the  replicate  test 
given  in  Reference  17,  i.e.,  Test  No.  1-10.  If  the  data  curve  in 
Figure  22  was  simply  shifted  to  the  right  a  small  amount,  the 
agreement  would  be  excellent.  Such  a  shift  appears  to  be  justified  on 
the  basis  of  the  data  scatter,  since  a  comparison  of  the  hydrostatic 
compression  curves  in  Figures  7  and  22  reveals  that  the  measured 
hydrostatic  response  in  Figure  22  is  stiffer  than  the  data  given  in 
Figure  7  to  which  the  modei  was  fit.  Aside  from  this  difference,  the 
qualitative  features  of  the  plain  concrete  behavior  are  obviously 
captured  very  well  by  the  model . 

5.2  SERIES  NO.  3. 

Test  Series  No.  3  was  designed  to  examine  the  behavior  of  plain 
concrete  to  circular  stress  paths  in  the  deviatoric  plane.  The  stress 
paths  consisted  of  monotonic  hydrostatic  loading,  then  proportional 
deviatoric  loading  along  the  triaxial  compression  path,  and  finally  a 
circular  load  path  in  the  deviatoric  plane.  In  the  present  study,  the 
load  paths  corresponding  to  Tests  3-3  and  3-4,  as  well  as  Tests  3-10, 
3-11  and  3-12  were  considered;  Tests  3-3  and  3-4  had  the  same  load 
path,  as  was  also  the  case  for  Tests  3-10  to  3-12. 

Figure  23  gives  the  responses  predicted  by  the  model  to  the 
circular  stressed  paths,  as  well  as  the  corresponding  data.  Again, 
the  overall  agreement  between  predictions  and  measurements  is 
considered  good  for  plain  concrete.  If  the  data  curve  in  Figure  23(b) 
were  translated  vertically  upward  by  a  sufficient  amount,  the  results 
would  show  very  good  agreement.  However,  such  a  translation  does  not 
appear  to  be  justified  if  one  considers  the  data  from  Tests  3-10  and 
3-11,  which  had  load  paths  identical  to  Test  3-12.  Thus,  the  model 
appears  to  be  too  soft  during  the  proportional  deviatoric  leg  of  Test 


The  results  shown  in  Figure  23  are  also  of  interest  from  another 
point  of  view.  Recently,  the  new  endochronic  theory  —  which  forms  the 
basis  of  the  model  developed  herein  —  was  criticized  by  Bazant, 
Krizek  and  Shieh^®^  as  being  unable  to  correctly  describe  the 
response  of  geomaterials,  such  as  concrete,  to  what  has  been  termed 
*  load) ng-to-the-s i de" .  Without  providing  supporting  experimental 
evidence,  they  claim  that  the  model  is  unrealistic  because  it  predicts 
purely  elastic  response  for  infinitesimal  load! ng-to-the-s ide.  To  the 
authors’  knowledge,  no  data  are  currently  available  for  plain  concrete 
which  define  its  response  to  I oadi ng-to-the-s i de.  Moreover,  while  the 
models’  response  to  infinitesimal  I oadi ng-to-the-s ide  is  purely 
elastic,  its  response  to  finite  I oadi ng-to-the-s ide  is  inelastic.  In 
Tests  3-3  and  3-12,  the  plain  concrete  was  subjected  to  loading-to- 
the-side,  at  least  at  the  beginning  of  the  circular  stress  path  and 
possibly  when  the  path  crossed  the  other  two  principal  stress  axes. 
In  view  of  the  reasonably  good  agreement  between  the  data  and  the 
predictions  for  these  paths,  the  criticisms  by  Bazant,  et  ai . (^9) 
appear  to  be  premature  and  unjustified.  For  further  elaboration  on 
this  point,  see  the  discussion  by  Read.^4) 


5.3  SERIES  NO.  4. 


This  series  consisted  of  the  largest  number  of  tests,  namely  22. 
In  general,  the  tests  were  designed  to  explore  several  aspects  of 
ax i symmetric  stress  paths,  i.e.,  paths  for  which  at  least  two  of  the 
principal  stresses  were  equal;  this  included  (a)  development  of  shear 
strain  due  to  hydrostatic  pressure  increments  in  the  presence  of  prior 
deviatoric  stress,  (b)  dependence  of  response  on  stress  path,  and 
(c)  behavior  under  nonproportional  deviatoric  stress  paths.  In  that 
which  follows,  eight  selected  stress  paths  from  Series  No.  4  which 


shed  light  on  each  of  the  above  aspects  of  plain  concrete  behavior  are 
considered. 

Test  4-1  and  its  replicate,  Test  4-2,  provide  insight  into  the 
development  of  shear  strain  during  hydrostatic  pressure  increments  in 
the  presence  of  non-vanishing  shear  stress,  i.e.,  shear-vo I umetr i c 
coupling.  These  tests  followed  the  stress  path  described  by  the 
insert  in  Figure  24.  This  figure  also  shows  the  response  to  the 
prescribed  stress  path  predicted  by  the  present  model,  as  well  as  the 
responses  measured  in  the  tests.  The  section  of  the  path  from  (b)  to 
(c)  is  of  primary  interest,  since  this  is  where  the  shear-volumetric 
coupling  occurs  under  purely  hydrostatic  increments.  As  the  results 
given  in  the  figure  indicate,  the  agreement  between  the  model  and  the 
two  sets  of  data  is  excellent. 

Next,  we  consider  Tests  4-9,  4-11,  4-13,  4-15,  4-16  and  4-17. 
This  group  of  tests  was  designed  to  explore  the  effect,  if  any,  of 
stress  path  dependence  on  plain  concrete  behavior.  Except  for  Test  4- 
13,  the  remaining  tests  consist  of  stress  paths  that,  while  quite 
different  from  one  another,  nevertheless  have  at  some  point  a  common 
stress  state,  namely,  =  8  ksi  and  rQ  =  4  ksi .  Thus,  by  comparing 
the  strains  measured  at  this  stress  state  in  the  various  tests  allows 
an  assessment  of  the  magnitude  of  the  stress  path  effect. 

Table  I  summarizes  the  principal  strains  in  the  various  tests 
that  were  measured  in  Reference  17  at  the  common  stress  state. 


Dependence  of  pressure  a  on  octahedra 


Table  1.  Summary  of  the  Principal  Strains  Measured 

at  the  Common  Stress  State  of  a  =8  ksi 

o 

and  x  =  ksi. 
o 


Test  Number 

Pr i nc i pa  1 

ex 

1  Strains  x 

ey 

103 

ez 

4-09 

5.10 

4.18 

15.89 

4-11 

2.86 

3.13 

11.15 

4-15 

1.64 

1.72 

12.72 

4-16 

10.25 

10.00 

13.27 

4-17 

1.19 

1.33 

12.78 

As  an  inspection  of  this  table  reveals,  there  is  substantial  depen¬ 
dence  of  the  deformation  at  the  common  stress  state  on  the  particular 
path  taken  to  get  to  that  state.  The  major  principal  strain,  ez, 

varies  from  11.15  x  10“3  to  15.89  x  10~3,  while  the  other  two  strains 

-3  -3 

show  much  greater  variation,  ranging  from  about  1  x  10”  to  10  x  10  . 

The  ability  of  the  present  model  to  describe  the  measured 
responses  for  Tests  4-9,  4-11,  4-13,  4-15,  4-16  and  4-17  is  illus¬ 

trated  in  Figures  25  to  32.  Unfortunately,  none  of  the  tests  in  this 
group  were  repeated,  so  we  have  no  insight  into  the  level  of  data 
scatter.  In  view  of  the  usual  scatter  found  in  concrete  data , however , 
the  agreement  between  the  model  predictions  and  the  data  is  considered 
excellent;  this  agreement  demonstrates  the  capability  of  the  present 
model  to  account  for  the  influence  of  the  stress  path  on  the  resulting 
response . 


PRINCIPAL  STRAINS  x  103 


Figure  27.  Measured  and  predicted  responses  for  Test  4-11. 
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The  final  ax i symmetric  test  to  be  considered  here  is  Test  4-20. 
The  prescribed  stress  path  for  this  test  consisted  of  nonproportional 
stress  deviation,  and  is  shown  in  Figure  33.  The  predicted  response 
to  this  prescribed  stress  path,  as  well  as  the  corresponding  data,  are 
also  shown  in  Figure  33.  Although  some  discrepancy  between  the  model 
and  the  data  is  evident  in  the  late  stages  of  the  stress  history,  the 
agreement  during  most  of  the  stress  history  is  excellent.  Therefore, 
the  overall  agreement  is  considered  good,  considering  the  usual  degree 
of  scatter  in  concrete  data;  this  test  was  also  not  repeated 
experiments  I ly . 

5.4  SERIES  NO.  5. 

The  purpose  of  the  Series  5  tests  was  to  investigate  the 
response  of  plain  concrete  to  unsymmetric  stress  paths  in  a  given 
deviatoric  plane.  In  the  present  study,  attention  was  focussed  on 
Tests  5-1  to  5-3  from  this  series,  with  emphasis  on  the  behavior 
during  the  stress  excursions  in  the  deviatoric  plane  corresponding  to 
aQ  s  4  ksi .  Details  of  the  prescribed  stress  paths  for  these  tests 
are  depicted  in  Figure  34(a).  The  measured  and  predicted  deviatoric 
strain  paths,  corresponding  to  those  portions  of  the  stress  paths  in 
the  deviatoric  plane  for  Tests  5-1  to  5-3,  are  shown  in  Figure  34. 
Also  included  in  this  figure  for  comparative  purposes  are  corres¬ 
ponding  predictions  made  by  Stankowsk i ,  using  the  Stankowsk i- 
Gerstle  model .  From  an  inspection  of  this  figure,  one  observes 
that  the  point  on  the  e2  axis  at  which  the  nonproportional  portion  of 
the  path  begins  for  the  present  model  is  consistently  higher  in  all 
three  tests  than  the  data.  It  appears,  therefore,  that  the  model  is 
too  soft  for  proportional  excursions  from  the  hydrostatic  axis  out 
along  the  tri axial  compression  axis,  confirming  the  same  observation 
made  earlier  in  conjunction  with  the  circular  stress  paths  in  Series 
3.  Despite  these  small  differences,  however,  Figure  34  shows  that  the 
present  model  does  an  excellent  job  of  predicting  the 
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Figure  34.  Stress  paths  for  Tests  5-1  to  5-3. 


measured  responses  in  the  deviatoric  plane,  and  is  clearly  superior  to 
the  Stankowsk i -Gerst I e  model . 

5.5  SERIES  NO.  6. 

The  tests  in  this  series  were  designed  to  investigate  the 
response  of  plain  concrete  to  piecewise  uniaxial  load  steps  which 
formed  closed  loops.  The  loading  was  applied  in  a  manner  such  that 
all  possible  uniaxial  and  equi-biaxiat  stress  states  of  a  given  level 
were  reached.  Attention  in  the  present  study  was  focussed  on  Test  6- 
2,  for  which  the  prescribed  stress  history  was  as  shown  in  Figure 
35(a).  The  peak  values  of  the  principal  stresses  for  this  path  were 

3.6  ksi . 

The  response  predicted  by  the  present  model  to  the  prescribed 
stress  path  is  illustrated  in  Figure  35(b),  where  the  corresponding 
data  are  also  given.  Considering  the  typical  scatter  in  concrete 
data,  the  agreement  between  the  model  prediction  and  the  data  is  again 
considered  excellent. 
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SECTION  6 


CONCLUSIONS 

An  advanced  nonlinear  constitutive  model  has  been  developed,  on 
the  basis  of  the  new  endochronic  plasticity  theory,  to  describe  the 
general  behavior  of  plain  concrete  over  the  range  of  stresses  where 
significant  cracking  does  not  occur.  The  model  is  isotropic,  rate- 
independent  and  thermodynamically  sound.  It  satisfies  Drucker’s 
Stability  Postulate,  Prager’s  Continuity  Condition  and  the  Second  Law. 
The  model  exhibits  shear-volumetric  coupling  in  the  form  of  shear 
compaction,  dependence  of  shear  behavior  on  hydrostatic  pressure, 
hysteresis  and  stress-path  dependence.  The  wide  variety  of  proof- 
tests  to  which  the  model  was  subjected  clearly  demonstrate  that,  over 
the  intended  range  of  application,  the  present  endochronic  model 
provides  an  excellent  description  of  the  behavior  of  plain  concrete. 

Our  future  research  activities  will  be  devoted  to  extending  the 
capabilities  of  the  present  model  to  a  wider  range  of  behavior;  this 
will  include  focus  on  the  following  tasks: 

•  Modify  the  unloading-reloading  shear  behavior 
through  the  addition  ox  shear  hardening  to  more 
closely  portray  experimental  data. 

•  Extend  model  to  cover  the  range  of  behavior  over 
which  significant  cracking  occurs. 


Introduce  strain-rate  dependence. 
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APPENDIX  A 


A  PROOF  OF  EQUATION  (3.50) 

The  purpose  of  this  appendix  is  to  show  that 


fz 


0 


where 


c  >  0,  M(z) 


p(z')dz' 

.  0 


and 


Um  =  M(«)  <  - 


Here,  p(z)  is  a  positive  monotonical ly  decreasing  function  of 
domain  0  <  z  <  *  and 


I  im 


/>(z) 


=  0 


Therefore,  M(z)  is  positive,  monoton i ca I  I y  increasing  and 
the  sense  that 


(A.l) 


(A.2a,b) 


(A.  3) 

z  i n  the 


(A. 4a) 

convex  in 


M(z)  -  M(z  -  a)  <  a  $ 


(A. 4b) 


Thus  in  vie*  of  Eqs .  (A . 2b) 

and  (A. 4) 

M(z)  -  M(z  -  a)  <  ap(z  -  a) 

(A. 4c) 

It  follows  from  Eqs.  (A. 4a) 

and  (A. 4c) 

that 

*iZ{M(z)  -  -  a)) 

=  0 

(A.  5) 

To  prove  the  validity  of 

theorem. 

Eq.  (A.l) 

we  first 

prove  the  following 

Theorem:  Let  R(z)  be  a 

positive  monoton i cal  ly 

decreasing  function 

bounded  from  above  and  below  in  the  sense  that 

0  <  R(z)  <  R(o) 

(A.  6) 

and 

^R(z)  =  0  . 

(A.  7) 

The  operational  definition 

of  Eq.  (A. 7) 

is  the  following:  Given  an  e 

however  small  there  exists 

a  z  <  *  such 

that 

R(z)  <  € 

(A. 8a) 

Then  the  following  is  true: 

z 

111!  />(z  -  z')R(z')dz* 

=  0 

(A. 8b) 

Jo 

IV 

tv 


.1 


rs 


The  integral  I  on  the  lefthand  side  of  Eq.  (A. 8b)  may  be  written  in 


the  form 


I  =  /J(z  -  z')R(z')dz'  ♦  p{ z  -  z')R(z')dz'  (A. 9) 

JO  Jz0 

where  0  <  z  <  z.  Because  of  the  properties  of  R(z),  it  follows  that 
o 


p{ z  -  z')R(z')dz'  <  R(o)  p(z  -  z')dz' 

0  Jo 


(A- 10) 


Or,  in  view  of  Eq.  (A. 2b), 


°p( z  -  z')R(z')dz'  <  R(o)  {m(z)  -  M[z  -  zq)} 

.  0 


(A- 11] 


Hence,  as  a  result  of  Eq.  (A. 5) 


p (z  -  z')R(z')dz'  =  0 


(A .  12] 


Also,  following  the  same  reasoning  as  above, 


p(z  -  z')R(z')dz*  <  R(z0)M(z  -  zo) 


(A. 13; 


APPENDIX  B 


NUMERICAL  SCHEME  FOR  TREATING  THE  PRESENT  CONSTITUTIVE 
MODEL  UNDER  STRESS-DRIVEN  CONDITIONS 


The  numerical  scheme  used  to  generate  the  model  predictions 
presented  in  this  report  is  described  below.  The  approach  consists  of 
an  explicit  incremental  Euler  scheme  developed  specifically  for  numer¬ 
ically  treating  the  governing  system  of  equations  for  the  model  under 
stress-driven  conditions.  Due  to  its  explicit  nature,  care  must  be 
taken  to  insure  that  the  prescribed  increments  are  sufficiently  small 
that  the  computed  behavior  is  independent  of  the  increment  size. 
Otherwise,  the  scheme  is  straightforward,  efficient  and  easy  to  imple¬ 
ment.  * 

The  basic  system  of  equations  which  govern  the  behavior  of  the 
model  are  as  fol lows: 


'(2s  *  2')dF  d2 


(B  1) 


*(2h  *  z'ld F  d2 


d2  =  2n (de  -  dgp) 


(B  2) 


(B  3) 


*.  This  numerical  scheme  is  due  to  Dr.  H.  Murakami,  a  consultant  to  S- 
CU8ED  and  Assistant  Professor,  University  of  California,  San  Diego. 
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6a  =  K (de  -  depJ 


(B .  4) 


dz2  =  dep  2  ♦  k2(dep) 2 


(B  5) 


dzs  =  dz/Fs 


(B .  6) 


dzH  =  dz/(kFH) 


(B .  7) 


Following  the  approach  described  in  Section  2,  the  kernel  functions 
p(zs)  and  1(zh)  are  represented  by  exponential  series  of  the  form: 


p  =  Hv 


(B8) 


*  *  ZbV''1” 


(B-9) 


It  then  follows  that  the  expressions  for  £  and  a  given  by  Eqs.  (B.l) 
and  (B.2)  can  be  alternately  written  as 


5  =  21  9r 


(B.10) 


=  x  p: 


(B  11) 


where  Qr  and  Pj  satisfy  the  following  ordinary  differential  equations: 


♦"  sT  Wi'  - 


From  Eqs.  (B.10)  to  (B.13),  we  can  write: 


dg  *  A  dep  -  9  dzg 
da  a  B  dep  -  P  dz^ 


where 


A  =  T]  A 
z—  r 

r 

S  =  I  *,  !r 

r 

B  =  E  Bi 

i 

p »  £  fii  pi 


(B ■ 14) 
(B  15) 


(B  16) 


Introducing  Eqs.  (B.6)  and  (B.7)  into  Eqs.  (B.14)  and  (B.15)  leads  to 
the  expressions: 


which  when  combined  with  Eq.  (B.5),  leads  to  the  following  quadratic 
expression  for  dz: 


2 

adz  ♦  b  dz  ♦  c  =  0 

where 


b 


I 


9  *  d5 


kP 


A2  F 


B2  Fl 


d  a 


d5  *  d5 

k  da 

2 

A2  * 

.  B  . 

(B  19) 


(B .  20) 


Thus,  for  prescribed  values  of  ds  and  da,  dz  is  given  by  the  equation: 


dz  s  *  2a  *  h  -  4  3C  (B  12) 

and  by  using  Eqs .  (B.17)  and  (B.18),  dep  and  dtp  can  then  be  found, 
after  which  de  and  de  are  determined  from  Eqs.  (B.3)  and  (B.4). 
Therefore,  this  approach  allows  the  strain  increments  d_e  to  be  found 
for  prescribed  increments  in  the  stress  do. 
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APPENDIX  C 


DERIVATION  OF  EQUATION  (4.12) 


When  a  i s  on  the  linear  portion  of  the  hydrostatic  stress-strain 
curve,  it  follows,  from  Eq.  (3.37),  that 


a  =  k  K 


d^H 

'SOdT 


(C.l) 


This  expression  may  be  combined  with  Eq .  (3.36),  i.e., 

dz2  =  df2  *  k2  d$^  (C . 2) 

to  give  the  following  result: 


dz^  1  - 


U1  *  * ) 


2)1/2 


=  d<T« 


(C  3) 


where  a ^  denotes  the  reference  pressure  on  the  linear  portion  of  the 
hydrostat . 

We  may  decompose  into  two  terms  and  write: 

Ko(‘  ♦  *«)  ■  K»(‘  *  ^h)  *  V?h  (C  4) 

where  is  the  value  of  <J||  due  solely  to  hydrostatic  compression, 
while  is  the  component  due  strictly  to  shear.  For  shear  at  a  =  c^, 
it  fol lows  that 


•i  *  ^C1  *  >  ?h) 


Hence,  Eq.  (C . 4)  can  be  rewritten  in  the  form: 


U1  ♦  p  Sh]  *  'l  *  V  ?H 


Upon  substituting  Eq.  (C.5)  into  Eq.  (C.3),  we  obtain  the  result 


=  1  - 


K  8 

1  *  ?H 


When  o  -  a^,  it  follows  from  Eq.  (C.l)  that 


»i  = »  ko(i  */>«■„•/>  ?;)tt 


so  that 


K  p  ) 

dz  =  k  1  *  £  d 

«.  1 


y  =  z  =  z 


K  p 

i  0  /-s 

x  *  1  ♦  77 
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where  z'  denotes  the  value  of  z  at  which  shearing  begins.  Then 
(C.8)  can  be  written 


dy  = 


ftojl 


k  a 

or 


x  dx 


Integrating  yields 

y 


_  1_  f  2 

"  2a  “ 


where  we  have  set 


From  Eqs.  (C . 6)  and  (C.10),  it  follows  that 


1/2 


(C 

(C 

(C 

(C 


Using  Eq.  (C.13),  we  can  rewrite  Eq.  (3.11)  in  the  form: 


dz  = 


x  dx 


(C 


since  dy  =  dz.  Eqs.  (C.14)  and  (C.15)  may  now  be  combined  to  give 


*r  =  1  («2  -  i)1'2  dx 
Ss  a 


(C 


which,  upon  integration,  leads  to  the  result: 


But  from  (C.12),  we  have 


x2  -  1  =  2  ay 
2 

x  =  1  ♦  2  ay 

which,  in  view  of  Eq.  (C.9),  may  be  rewritten  in  the  form: 
x2  -  1  s  2a  2s 
x2  =  1  ♦  2a  zs 
since  for  F^  =  1,  we  have 
dzs  *  dz 


(C.18) 


(C.19) 


and 


=  z  - 


z' 


(C.20) 


Upon  combining  Eqs.  (C.17)  and  (C.19),.  the  following  expression  is 
obta i ned : 


1  ♦  2a  z 

s 


log  2a  zg  ♦  1  ♦  2a  zg  j j  (C.21) 


which  is  Eq.  (4.12)  given  in  the  main  text. 
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